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Abstract 

Highly accurate results are presented for the susceptibility, x(T) of the s = 
1/2 Heisenberg antiferromagnetic chain for all temperatures, using the Bethe 
ansatz and field theory methods. After going through a rounded peak, x{T) 
approaches its asympotic zero-temperature value with infinite slope. 
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In a pioneering work in 1964, Bonner and Fisher estimated numerically^ the susceptibility, 
x(T) for the s = 1/2 Heisenberg antiferromagnetic chain, using chain lengths of up to 11. 
The exact value at T = was obtained using the Bethe ansatzi'i. Since that time, the 
Bonner-Fisher curve has frequently been used by experimentalists to determine the value of 
the exchange coupling, J, and to determine whether or not a given material has sufficiently 
small anisotropic exchange and inter-chain couplings to be approximated by this model. 
Recently it has become possible to calculate this curve much more accurately, using the 
Bethe ansatzE This method easily gives very accurate results for T not too small, but 
becomes increasingly difficult as T — > 0. On the other hand, an analytic formula can be 
derived for x{T) at small T, from conformal field theory methods. Perhaps surprisingly, 
this formula predicts, due to the marginally irrelevant operator, that the susceptibility has 
infinite slope at T = 0. Here we present results for x(T) obtained using the Bethe ansatz for 
T > .003J and compare them to the conformal field theory prediction, obtaining excellent 
agreement. The field theory prediction for the general xxz model is also given. Figure (1) 
shows the susceptibility obtained from the Bethe ansatz and Figure (2) compares this result 
to the field theory prediction at low T. Note that, with decreasing T, after passing through 
the maximum, \ ~ -147/ J, at T « .640824J, the slope of \ starts to increase below the 
inflection point at T « .087 J, approaching oo as T — > 0. 

The Heisenberg Hamiltonian is written: 

H = JJ2S. t -S, l+1 . (1) 

i 

The susceptibility per site is given by: 

xcn^E<^>, (2) 

i 

The low energy effective field theory description! is given by the k — 1 Wess-Zumino-Witten 
(WZW) non-linear cr-model, or equivalently a free, massless boson, with an effective "velocity 
of light" or spin-wave velocity: 

v = txJ/2. (3) 

This value of v is determined from the slope of the dispersion relation obtained from the 
Bethe ansatz. The uniform part of the spin density is given by the conserved current 
operators, Jl, Jr for left and right-movers: 

Si w J L {xi) + J R {xi). (4) 

In the WZW model Jl and Jr are uncorrelated and their self-correlations are given by: 
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This result can be extended to finite temperature by a conformal transformation. This 
simply replaces: 



v(3 

v r =F ix — ► — sin 
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where (3 = 1/T. 

The susceptibility in the WZW model is thus given by: 
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The integral can be done by the change of variables: u = tan ^ and w = — 2 tan giving: 

(8) 
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Note that the integral is independent of r or u as it should be since the total spin is conserved. 
Plugging in the spin- wave velocity, v = Jn/2, gives the zero-temperature susceptibility, 
\(0! 1 ./r; 2 . 

The fact that x{T) is independent of T in the WZW model follows from scale invariance. 
To obtain the leading T-dependence, we must perturb about the scale invariant fixed point 
Hamiltonian with the leading irrelevant operator. This perturbation is written: 



M , -8tt 2 v 7 7 
on = — —gJL ■ Jr- 



v/3 



(9) 



This term is marginally irrelevant for the coupling constant, g > 0. We wish to calculate the 
correction to x(T) to first order in b~7i. Expanding e~ J d,xdT5H(T,x) ^ Q g rg ^ orc [ er ^ e correction 
to x involves four current operators. Due to the fact that the left and right currents are 
uncorrelated, this expression factorizes into a product of two two-point Green's functions, 
one for left-movers and one for right-movers. Using translational invariance, the spatial 
integrals factorize into two independent integrals of the form of Eq. (pi) , giving: 
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Again, the correction is naively temperature- independent, since g is dimensionless. However, 
this formula can be improved by replacing g by g(T), the effective renormalized coupling at 
temperature T. By integrating the lowest order (3 -function, g(T) is given by:! 



g(T) 
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47r<7iln(Ti/T)/V3' 



(11) 



Here g\ is the value of the effective coupling at some temperature T±. Both g\ and g(T) 
must be small for this formula to be valid. We may write this more compactly, as: 



g(T) 



4vrln(To/T) : 



(12) 
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for some temperature To. Note that the asympotic behaviour at small T becomes indepen- 
dent of <7i or T up to a correction of 0((lnT)~ 2 ). We expect higher order corrections to 
x(T) of 0(g(T) 2 ), 0(g(T) 3 ), etc. By shifting T , we can eliminate the O(lnT) -2 ) part of 
the 0(g(T) 2 ) correction. Thus we obtain the leading T-dependence of x- 

^^ + S?tok7r) +0((1,irr3) ' (13) 

This formula is universal in the sense that if we add additional interactions to the Heisenberg 
Hamiltonian that respect the SU(2) and translational symmetry and are not sufficiently 
large as to drive it into another phase, then this formula applies for some values of v and 
To. In particular, if we add an antiferromagnetic second nearest neighbour interaction, J2, 
the value of the bare coupling, g, decreases and reaches zero at J2 ~ .24J. At this point all 
logarithmic terms in x{T) vanish, corresponding to T — > 00 and x{T) should be analytic 
near T = 0. For larger J 2 the system develops a gap and x{T) vanishes exponentially as 
T — > 0. Eq. (|13"D should be valid for arbitrary half-integer spin Heisenberg antiferromagnets. 
For the ordinary s — 1/2 Heisenberg model, v = Jtt/2 giving: 

J ^ T)xl + 2^km- (14) 

A good fit to the Bethe ansatz data is obtained with To ~ 7.7 J, as shown in Figure (2). 
Eq. (|TJ]) is valid to within 2% for T < .1J. Similar formulas for the finite-size dependent 
of low energy states were obtained in Ref. ( ||) with T /T replaced by L/L where L is 
the size of the system. The role of infrared cut-off on the renormalization of the coupling 
constant is played by either the length, L or an effective thermal length v/T in the two cases. 
Note that the field theory calculation of the susceptibility is done in the limit LT/v — > 00 
whereas the finite-size spectrum calculations are done in the opposite limit LT/v — * 0. In 
both cases the space-imaginary time surface is an infinite length cylinder of circumference 
v/T or L respectively. By pushing all calculations to second order in g, predictions could be 
made relating the various values of the Lo's for different energy levels and Tq. Alternatively, 
we may use the deviation of x{T) and the energy levels from their asympototic values to 
define four different estimates of the effective coupling constant versus length. These are 
shown in Figure 3. The singlet and triplet excited state estimates of the effective coupling, 
g(L), only converge very slowly as g{L) 2 , as we expect since, in general, all these quantities 
receive corrections at next order in g. Surprisingly, the susceptibility and triplet estimates 
of g, using the relationship L <-> v/T, appear to converge much more rapidly as 1/L 2 , 
suggesting the absence of corrections to any finite order in g(L). [The groundstate and 
triplet estimates appear to converge rapidly to a small non-zero difference. This remains a 
puzzlingly descrepancy between conformal field theory and Bethe ansatz calculations. See 
Ref. ( 0) for further discussion.] 

We now consider the general xxz model: 



# = J £ [S?S* +1 + SfSf +1 + AS*S? +1 ] . (15) 



For A > 1, the Hamiltonian has a Neel ordered groundstate and a gap. Hence the sus- 
ceptibility vanishes exponentially at low T. For A < 1, the system remains gap less. It is 
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now convenient to use abelian bosonization, involving a free boson 0. The zz component 
becomes: 



1 

47T 



130 

v dt 



dx 



The other components of the interaction become: 

JlJr + J V lJr<x cos V8^(f>. 
The zz part is exactly marginal and has the effect of rescaling the boson:i 



f^KR 



with R < 1/ v27r. The other part then becomes: 



J£J£ + «occos(20/ J R). 



This has scaling dimension 



x 



1/tvR 2 > 2 
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and is irrelevant. It is the leading irrelevant operator, provided x < 4. A renormalization of 
the velocity also occurs. The zero temperature susceptibility gets rescaled to: 



x(o) 



v(2nRy 



(21) 



The first order contribution to the susceptibility from cos(20/i?) vanishes. The second order 
contribution is determined by a standard scaling argument, giving: 



X(T) 



v(2nR) 2 



+ constant -T^ 2 " 4 ). 



(22) 



Again this formula is universal in the sense that the two parameters v and R determine all 
low energy features of the model and R < 1/ \^2n. Increasing anisotropy leads to decreasing 
R. Eq. ( |2"2] ) should apply for arbitrary half-integer spin xxz antiferromagnets in the gapless 
phase, provided that 1/nR 2 < 3. Otherwise the exponent is replaced by 2. For the s = 1/2 
Heisenberg model -R(A) and v (A) have been determined from the Bethe ansatz.Hi'1 Letting: 



A = cos 9 



(23) 



we find: 



2nR(A) = yi- - 

Jtt sin 6 



v(A) 



20 



(24) 
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Thus 



X(0) = (25) 

7r(7r — V) sm& 



(in agreement with the explicit Bethe ansatz calculations) and the exponent is given by 



2/^-4 = -^-. (26) 

7T — t7 



For small anisotropy: 



, 2 , 4 V 2 (1 " A) 



2tt/R 2 - 4 « -. (27) 



7T 



x(T) has an infinite slope at T = for A > cos 1 (7r/5) ~ .809. For A < .5, the exponent 
(2n/R 2 - 4) is replaced by 2. 

These results may prove useful in experimental studies of quasi-one-dimensional anti- 
ferromagnets. In particular, we may define a "critical region" where the temperature de- 
pendence is governed by the leading irrelevant operator so that Eqs. (|13|) or (|22|) hold. 



For the isotropic Heisenberg model, this critical region can be identified approximately as 
T < .1 J, below the inflection point. Only if the Neel temperature, determined by interchain 
couplings, is below this value can the one dimensional critical behaviour be observed. 
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FIGURES 



FIG. 1. x(T) from the Bethe ansatz. x(0) is taken from Ref. ( ||). 
FIG. 2. Field theory [Eq. (pi)] versus Bethe ansatz formulas for x(T) at low temperature. 



FIG. 3. Estimates of the effective coupling from lowest order perturbation theory correction 
to finite-size energy of groundstate, first excited triplet state, first excited singlet state and to 



susceptibility, using L *-* v/T. The renormalization group prediction of Eq. (12) is also shown. 
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